Effects of confinement on permanent electric-dipole moment of Xe atoms in liquid Xe 
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Searches for permanent electric-dipole moments (EDM) of atoms provide important constraints on 
competing extensions to the standard model of elementary particles. Recently proposed experiment 
with liquid 129 Xe [M.V. Romalis and M.P. Ledbetter, Phys. Rev. Lett. 87, 067601 (2001)] may 
significantly improve present limits on the EDMs. To interpret experimental data in terms of 
CP-violating sources, one must relate measured atomic EDM to various model interactions via 
electronic-structure calculations. Here we study density dependence of atomic EDMs. The analysis 
is carried out in the framework of the cell model of the liquid coupled with relativistic atomic- 
structure calculations. We find that compared to an isolated atom, the EDM of an atom of liquid 
Xe is suppressed by about 40%. 
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Most extensions of the standard model of elementary 
particles, e.g., supersymmetry, naturally produce per- 
manent electric dipole moments (EDM) of atoms and 
molecules that are comparable to or larger than 
present limits (see, e.g., a popular review |2|). For 
example, the most accurate to date determination of 
atomic EDM of 199 Hg Q sets limits on a number of im- 
portant parameters: CP-violating QCD vacuum angle, 
quark chromo-EDMs, semileptonic CP-violating param- 
eters, and restricts parameter space for certain extensions 
to the standard model. A substantial, several orders of 
magnitude improvement in sensitivity to all the enumer- 
ated sources of CP-violation is anticipated in an exper- 
iment proposed by Romalis and Ledbetter [j]. These 
authors propose to search for an EDM of a liquid sample 
of 129 Xe. Compared to the gas-phase experiment |5j, a 
drastically improved sensitivity of the liquid Xe experi- 
ment is mainly due to the higher number densities of the 
liquid phase (10 22 atoms/cm 3 ). 

The very use of the liquid phase raises questions about 
density-dependent factors which can influence the out- 
come and interpretation of the experiment. For exam- 
ple, an EDM experiment with a molecular liquid was 
proposed in Ref. The authors found an additional 

suppression of the EDM signal by a factor of a hundred 
due to a reduced population of molecular rotational levels 
in liquid. Although the experiment with liquid Xe will be 
free from such an effect, it is clear that the effects of the 
liquid phase on atomic EDMs have to be investigated. 

An EDM of an atom is related to a strength of a CP- 
violating source via electronic-structure (enhancement or 
shielding) factors. For an isolated Xe atom such factors 
were computed previously: P,T-odd semileptonic inter- 
actions were considered by Martensson-Pendrill and 
nuclear Schiff moment — by Dzuba et al. 8]. Here we 
employ a simple cell model to study density dependence 
of the electronic-structure factors. Technically, we ex- 
tend the previous atomic relativistic many-body calcula- 
tions by confining a Xe atom to a spherically-symmetric 
cavity. In a non-polar liquid such as liquid Xe, this cav- 
ity roughly approximates an averaged interaction with 



the neighboring atoms. Imposing proper boundary con- 
ditions at the cavity radius, first we solve the Dirac- 
Hartree-Fock (DHF) equations and then employ the rela- 
tivistic random-phase approximation (RRPA) to account 
for correlations. To the best of our knowledge, here we 
report the first ab initio relativistic calculations of prop- 
erties of a liquid. We find that compared to the EDM of 
an isolated atom, the resulting EDM of an atom of liquid 
Xe is suppressed by about 40%. Thus if the experiment 
with liquid Xe is carried out with the anticipated sensi- 
tivity, we expect that the inferred constraints on possible 
sources of CP-violation would be indeed several orders of 
magnitude better than the present limits. 

Sources of atomic EDM — The conventional atomic 
Hamiltonian Hq among other symmetries is invariant 
with respect to space- reflection (P) and time reversal 
(T). Therefore, on very general grounds, an expectation 
value of the electric dipole operator D — — ^ in a 
non-degenerate atomic state \ ^o) vanishes. The tiny CP- 
violating interactions, here generically denoted as Hqp = 
J"V /icp( r i)i break the symmetry of the atom and induce 
a correction to the electronic state \S&) = |^o) + |^) ■ To 
the lowest order 

i**)=I>*> Egggga . a) 

k 

where Ek and \^k) are eigenvalues and eigenfunctions 
of Ho. Due to selection rules, the \6^) admixture has 
a parity opposite to the one of the reference state I'J'o). 
Because of this opposite-parity admixture the atom ac- 
quires a permanent EDM 

d= (*|D|*) =2(tf |D|J*). (2) 

Now we specify particular forms of Hep- An anal- 
ysis pj shows that for diamagnetic atoms, such as Xe, 
the EDM predominantly arises due to P,T-odd semi- 
leptonic interaction i?TN between electrons and nuclc- 
ons and also due to interaction i?sM_of electrons with 
the so-called nuclear Schiff moment 9] . Smaller atomic 
EDM is generated by intrinsic EDM of electrons and 
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we will not consider this mechanism here. Atomic units 
|e| = h = m e — 47T£ = 1 are used throughout. 

Explicitly, the effective P,T-odd semilcptonic interac- 
tion Hamiltonian may be represented as [jj 



/*TN(r e ) = V2Gf Ctn cr N ■ (17075 cr) e pn (r e ) 



(3) 



Here subscripts e and N distinguish between operators 
acting in the space of electronic and nuclear coordinates 
respectively. Ctn is a coupling constant to be determined 
from an interpretation of EDM measurements and to be 
compared with theoretical model-dependent predictions. 
Due to averaging over nuclear degrees of freedom, this in- 
teraction depends on nuclear density distribution Pn{t)- 
In the following, we approximate pn{t) as a Fermi distri- 
bution pn{t) = Po/(l + cxp[(r — c)/a]) with c = 5.6315 
fm and a = 0.52 fm. Finally, G F w 2.22254 x 1CT 14 a.u. 
is the Fermi constant. 

The interaction of an electron with the nuclear Schiff 
moment S has the form [Toj 



hsM{r e ) = -tt Pn (r e ) (r e • S) , 

i>4 



(4) 



where B4 — J Q °° r^pM (r) dr is the fourth-order moment 
of the nuclear distribution. The Schiff moment character- 
izes a difference between charge and EDM distributions 
inside the nucleus. It depends on a number of important 
CP-violating parameters enumerated in the introduction. 

Finally, we emphasize that both /itn and /ism are con- 
tact interactions. They occur when an electron pene- 
trates the nucleus. The electron speed at the nucleus is 
approximately aZc ~ \c (Z = 54), i.e., a fully relativis- 
tic description of electronic motion is important in this 
problem. 

Cell model of liquid xenon — Here we employ a simple 
cell model (see [llj and references therein) to estimate 
the effects of the environment on permanent EDM of a 
given atom. According to the cell model, we confine an 
atom to a spherical cavity of radius 



p i 3 1 

47T n 



1/3 



(5) 



n being the number density of the sample. For a den- 
sity of liquid Xe of 500 amagat [T^. R ca , v — 4.9 bohr. 
In non-relativistic calculations periodicity requires that 
the normal component of the gradient of electronic wave- 
function vanishes at the surface of the cell (see, e.g., 01) 



dr 



0. 



(6) 



Before proceeding with a technical question of imple- 
menting these boundary conditions in relativistic calcu- 
lations, we notice that the cell model implicitly incor- 
porates an average polarization interaction with the me- 
dia. Indeed, the Hamiltonian of an atom placed in the 
liquid in addition to the conventional atomic Hamilto- 
nian Hq includes interaction of electrons with the rest of 



the atoms in the media. This interaction is dominated 
by polarization potential. An important point is that 
the averaged polarization interaction can be expressed as 
V p = —1/2(1 — e _1 )i? c T a ^, where e is the dielectric con- 
stant of the media |13| . This interaction does not depend 
on electronic coordinate — it is just an additive constant 
which does not affect calculations of EDM. Thus we may 
approximate the total Hamiltonian with the traditional 
atomic Hamiltonian Hq. 

Further, the spherical symmetry of the cell allows us 
to employ traditional methods of atomic structure. The 
only modification is due to boundary conditions 10. 
However, in relativistic calculations, a special care should 
be taken when implementing this boundary condition. 
Indeed, the Dirac bi-spinor may be represented as 



(7) 



where P and Q are the large and small radial components 
respectively and f2 is the spherical spinor. The angular 
quantum number k = {I — j) (2j + 1). The nonrclativis- 
tic boundary condition © applied directly to the above 
ansatz would lead to two separate constraints on P and 
Q. This over-specifies boundary conditions and leads to 
the Klein paradox. 

A possible relativistic generalization of the boundary 
condition (JSJ) is 



d_l\ 

dr r 



■{Rc 



d QnK 
dr r 



{Rcav) ■ 



(8) 



Since in the non-relativistic limit the small component 
Q vanishes, this generalization subsumes Eq. @. Due 
to a semi-qualitative nature of our calculations, here we 
have chosen to use simpler (MIT bag model) boundary 
condition 



P71K {Rc 



Qtlk {Rc 



(9) 



Non-relativistically it corresponds to impenetrable cav- 
ity surface. Compared to this condition, the periodic 
boundary conditions JSJl are "softer", i.e., they modify 
the free-atom wavefunctions less significantly; we expect 
that our use of Eq. © would somewhat overestimate the 
effects of confinement in the liquid. 

Atom in a cavity: DHF and RRPA solutions — To re- 
iterate the discussion so far, within the cell model, the 
complex liquid-structure problem is reduced to solving 
atomic many-body Dirac equation with boundary condi- 
tions ©. The atomic-structure analysis is simplified by 
the fact that Xe is a closed-shell atom. Below we self- 
consistently solve the DHF equations inside the cavity. 
Then we employ more sophisticated RRPA. 

At the DHF level, the atomic wavefunction is repre- 
sented by the Slater determinant composed of occupied 
(core) orbitals tp a . These orbitals are determined from a 
set of DHF equations 



(c(a ■ p) + /3c 2 + V nuc + Vdhf) fa 



£afa 



(10) 
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where V nuc is a potential of the Coulomb interaction with 
a finite-size nucleus of charge density (r) and Vdhf is 
non-local self-consistent DHF potential. The DHF poten- 
tial depends on all the core orbitals. Similar equations 
may be written for (virtual) excited orbitals ip m . 

We solved the DHF equations in the cavity using a B- 
spline basis set technique by Johnson et al. [ljj. This 
technique is based on the Galerkin method: the DHF 
equations are expressed in terms of an extremum of an 
action integral Sa- The boundary conditions are incor- 
porated in the Sa as well. Further, the action integral 
is expanded in terms of a finite set of basis functions 
(B-splines). Minimization of such Sa with respect to 
expansion coefficients reduces solving integro-differential 
DHF equations to solving symmetric generalized eigen- 
value problem of linear algebra. The resulting set of basis 
functions is finite and can be considered as numerically 
complete. In a typical calculation we used a set of basis 
functions expanded over 100 B-splines. 

Given a numerically complete set of DHF eigenfunc- 
tions {<fi}, the permanent atomic EDM, Eq.Q, may be 
expressed as 



jDHF 



{<Pa\r\<Pm){<Pm\hcp\<Pa) 



(11) 



where a runs over occupied and m over virtual orbitals. 
Here hep is either a semileptonic interaction, Eq. J3J), or 
an interaction with the nuclear Schiff moment, Eq. 
An additional peculiarity related to the Dirac equation is 
an appearance of negative energy states (e m < — m e c 2 ) 
in the summation over intermediate states in Eq. 1)11)1. 
We have verified that these states introduce a completely 
negligible correction to the computed EDMs. 

To improve upon the DHF approximation, we have also 
computed EDMs using RRPA method 0] . This approx- 
imation describes a dynamic linear response of an atom to 
a perturbing one-body interaction (e.g., Hqp). The per- 
turbation modifies core orbitals thus changing the DHF 
potential. This modification of Vdhf m turn requires the 
orbitals to adjust self-consistently. Such a readjustment 
process defines an infinite series of many-body diagrams, 
shown, e.g., in Ref. The RRPA series can be summed 
to all orders using iterative techniques or solving DHF- 
like equations. We used an alternative method of solu- 
tions based on the use of basis functions ^(| . As an input, 
we used the DHF basis functions generated in the cavity 
(see discussion above), i.e. the boundary conditions were 
satisfied automatically. As a result of solving the RRPA 
equations, we have determined a quasi-complete set of 
particle-hole excited states and their energies. Then the 
EDMs are determined using expressions similar to Eq. 
and ©. 

Discussion and conclusions — First, we present the 
results of our calculations for an isolated atom (i? cav = 
oo). For the Schiff-moment-induced EDM, our results, 
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are in agreement with the recent calculations by Dzuba 
et al. @|. For the EDM induced by T,P-odd semileptonic 
interactions we obtain 



jDHF 
2 TN 



= 8.44 x 10 



d™ PA = 10.7 x 10" 



-13 r> 
-13 r 



TN &N a.u. 
TN fjv a.u. 



These values are to be compared with the results by 
Martensson-Pendrill 7|, d^ F = 7.764 and d^ PA = 
9.808 in the same units. The reason for the 10% dif- 
ference between our results and those from Ref. |7J is not 
clear. 

Before presenting results for finite cavity radii, let us 
consider individual contributions to EDM from various 
shells of Xe atom. These contributions for the Schiff- 
moment-induced EDM of an isolated atom are listed 
in Table Q] A similar table, but for the EDM arising 
from semileptonic interactions is given in Ref. 0. From 
these tables we observe that the dominant contribution to 
EDMs comes from the outer n = 5 shell. Thus we antic- 
ipate that a noticeable density dependence should occur 
when i? cav becomes comparable to the size of external 
n = 5 shell. We also notice that the contribution from 
the outer shell is relatively more important in RRPA cal- 
culations than at the DHF level, i.e., the RRPA results 
should exhibit stronger density dependence. 





DHF 


RRPA 


n = 1 


0.039 


0.039 


n = 2 


0.091 


0.092 


n — 3 


0.20 


0.21 


n = 4 


0.52 


0.64 


n = 5 


2.0 


2.8 


Total 


2.88 


3.78 



TABLE I: Individual contributions from various shells to the 
EDM of a free 129 Xe atom in the DHF and RRPA methods. 
The EDM is induced by the nuclear Schiff moment and it is 
given in units of S/(efm 3 ) x 10 _ls ecm. 



These qualitative conclusions for a confined atom are 
supported by our numerical results, presented in Fig. ^ 
Here we plot the ratios of atomic EDMs for the confined 
and isolated atoms as a function of i? C av ■ The EDMs be- 
come smaller as the density increases, n cx R~^- At the 
density of liquid Xe, i? cav ~ 4.9 bohr , the more accu- 
rate RRPA results show a 25% suppression of the atomic 
EDM due to confinement. Overall there is a noticeable 
density-dependence of atomic EDM. We expect the EDM 
signal (if found) to be broadened. The relevant charac- 
teristic width of the signal can be simply estimated from 
our Fig. ^ from the mean density fluctuations. 

^From Fig. ^ we notice that both semileptonic- and 
Schiff-moment-induced EDMs scale with i? cav in a sim- 
ilar fashion. This similarity can be explained from the 
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following arguments. The values of CP-violating matrix 
elements, Eq.© and Eq.QJ, are accumulated inside the 
nucleus. Non-relativistically, as r — > the wavefunc- 
tions scale as <p n i m {r) « N n i(R cav ) x r l Yi m (r), where 
N n i are normalization factors. Therefore the dominant 
contribution to the EDM, Ea. (|ll|) arises from mixing of 
s and p states. By factorizing the matrix element of hep 
as ((pns\h C p\<Pn'p) ~ N ns (R ca , v )N n ' p (R ca . v ) x (s\h C p\p) 
we see that the i? cav -independent factor (s\hcp\p) can 
be pulled out of the summation over atomic orbitals in 
Ea. (|ll(l . Thus, both semileptonic- and Schiff-moment- 
induced EDMs exhibit approximately the same scaling 
with the cavity radius. A correction to this "similarity 
scaling law" may arise, for example, due to different se- 
lection rules involved for the two EDM operators. 



It is worth emphasizing the semi-qualitative nature 
of our calculations. The analysis can be improved by 
employing more realistic models of liquid environment. 
Even within the cell model we could further refine our 
analysis. A dense liquid may be considered as a solid 
with vacancies, i.e., the decrease of the average bond 
length is negligible, rather the nearest-neighbor occu- 
pation numbers are decreased compared to solid. Xe 
condenses into face-centered cubic structure. The first 
nearest-neighbor shell contains twelve atoms (partially 
justifying the spherical symmetry of the elementary cell). 
The density of the solid Xe is 3.54 g/cm 3 , implying the 
half-radius of this shell of 4.2 bohr, somewhat smaller 
than i? C av ~ 4.9 bohr for liquid Xe. As follows from 
Fig. ^ this difference leads to more pronounced suppres- 
sion of the atomic EDM by 40% per cent. 



FIG. 1: The ratios of atomic EDMs for the confined and iso- 
lated atoms (suppression factor) as a function of cavity radius. 
The upper and lower sets of two curves are obtained with 
the DHF and RRPA methods respectively. EDMs induced 
by P,T-odd semileptonic interactions are shown as solid and 
dashed lines, while EDMs due to the Schiff moment — as 
dotted and dashed-dotted lines. Heavy dot marks our final 
results for liquid Xe. 



To reiterate, our work was motivated by anticipated 
significant improvements in sensitivity to atomic EDMs 
in experiments with liquid 129 Xe |4|. Here we investi- 
gated confining effects of the environment on the EDM 
of Xe atom. We carried out the analysis in the frame- 
work of the cell model coupled with relativistic atomic- 
structure calculations. We found that compared to an 
isolated atom, the EDM of an atom of liquid Xe is re- 
duced by about 40%. Thus if the experiment with liquid 
Xe is carried out with the anticipated sensitivity, we ex- 
pect that the inferred constraints on possible sources of 
CP-violation would be indeed several orders of magni- 
tude better than the present limits. 
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